Program for calculating bounds on the minimum rank of a graph 

using Sage 

Laura DeLoss* Jason Grout* Tracy McKay* Jason Smith* Geoff Tims* 

December 9, 2008 

Abstract 

The minimum rank of a simple graph G is defined to be the smallest possible rank over all 
symmetric real matrices whose ijth entry (for i ^ j) is nonzero whenever is an edge in 

G and is zero otherwise. Minimum rank is a diflicult parameter to compute. However, there 
are now a number of known reduction techniques and bounds that can be programmed on a 
computer; we have developed a program using the open-source mathematics software Sage to 
implement several techniques. In this note, we provide the source code for this program. 

Keywords, minimum rank, maximum nullity, zero forcing number. Sage program, mathematical 
software, symmetric matrix, rank, matrix, tree, planar graph, graph. 
AMS subject classifications. 05C50, 15A03 

1 Introduction 

In this note, we provide a listing for our Sage O version 3.1.2] program which computes upper and 
lower bounds for the minimum rank of a graph. 

We hrst include several example Sage sessions, illustrating the main function of the program, 
minrank_bounds. The lines starting with sage: are input lines (the "sage:" is not typed), while 
the other lines are output lines. Without any options, the niinrank_bounds function returns two 
numbers: a lower bound and an upper bound for the minimum rank. In the following example, we 
find that K'^ has minimum rank 1 (i.e., the lower bound and the upper bound are both 1). 

sage : minrank_bounds (graphs . CompleteGraph (3) ) 
(1, 1) 
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With the all_bounds option set to True, the function returns two hstings in the format "name: 
value" , where the name is the name of a bound and the value is the upper or lower bound. In the 
following example, we calculate bounds on the minimum rank of the Petersen graph. In the second 
call, we see that lower bounds are minimally given by the obvious bound on ranks (i.e., the rank 
has to be at least zero) and also by a minimal zero forcing set (a lower bound of 5). For upper 
bounds, we have a clique cover which gives a bound of 15, the fact that the graph is not outer 
planar (upper bound of 7), the fact that the graph is not a path (upper bound of 8), the fact that 
the graph is not planar (upper bound of 6), the fact that the minimum rank is at most |G| — 1 
(upper bound of 9), and the trivial upper bound on ranks (i.e., the rank is at most |G|). 

sage : minrank_bounds (graphs . PetersenGraphO ) 
(5, 6) 

sage: minrank_bounds (graphs. PetersenGraphO , all_bounds=True) 

(-['rank': 0, 'zero forcing': 5}, 

{'clique cover': 15, 

'not outer planar': 7, 

'not path' : 8, 

'not planar ' : 6 , 

'order' : 9, 

'rank': 10» 

If the tests option is set to a list, then only those tests are run. In the following example, we 
compute bounds on the minimum rank of the Heawood graph by using the zero forcing test and by 
testing if the graph is not planar. The "rank" tests are always run (and always give for a lower 
bound and \G\ for an upper bound). 

sage: minrank_bounds (graphs . HeawoodGraphO , all_bounds=True , 
sage: tests=['zero forcing', 'not planar']) 

({'rank': 0, 'zero forcing': 8}, {'not planar': 10, 'rank': 14}) 

For more complete documentation and a list of tests that can be run, as well as for sev- 
eral more examples, print the help by typing the function name followed by a question mark: 
minr aiik_bounds ? . 

This program contains and uses the minimum rank data listed in [2]. 

2 Program listing 

In the typeset program listing below, some lines are automatically broken that are not actually 
broken in the source code. If a line is broken into two lines in the listing, but should appear as one 
line in the program, then the line will end with and the remainder of the line will start with 
To use this program, download and extract the "source" from arxiv.org. Then you can 
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• upload the accompanying minrank . sws Sage worksheet into a notebook, or 

• copy the contents of the accompanying minrank. sage file into a cell of a Sage worksheet, or 

• load the minrank. sage file into a running Sage terminal session. 
We now proceed with the program listing. 

# Imports all the graphs of order 7 or less and stores # 

# them in a list called atlas.graphs so that # 

# atlas_graphs [ i ] is the ith graph in the atlas of # 

5 # graphs # 

/////////////////////////////////////////////////////////////////////////////////////^^ 



import networkx . generators . atlas 
10 atlas_graphs = [Graph(i) for i in \ 

networkx . generators . atlas . graph_atlas_g () ] 

#A list "database" of all the minimum ranks of graphs # 
15 # of order 7 or less . 



# 

# The minimum ranks are stored as ordered pairs : the 

# first coordinate is the graph number, the second 

# coordinate is the minimum rank. The first tuple in 

# the list min.ranks is just a position holder. 

/////////////////////////m^^^^ 



# 
# 
# 

# 
# 
# 



min_ranks 



(7,1) 
(16,2 
(24 
(32 
(40 
(48 
(56 
(64 
(72 
(80 



(96 



2) 
2) 
3) 
2) 
2) 
2) 
3) 
4) 
3) 
4) 



(8,0) 
(17 
(25 
(33 
(41 
(49 
(57 
(65 
(73 
(81 
(89 
(97 



[(0,None), ( 
, (9,1) , (10 

2) , (18,1) 

3) , (26,3) 

2) , (34,3) 

3) , (42,2) 

2) , (50,2) 

1) , (58,2) 

3) , (66,4) 

3) , (74,3) 

4) , (82,3) 

2) , (90,3) 
4) , (98,4) 



1 ,0) , (2 
,2), (11 
(19,0) 
(27,2) 
(35,3) 
(43,3) 
(51 ,2) 
(59,3) 
(67,2) 
(75,3) 
(83,5) 
(91,2) 
(99,4) 



,0) , (3, 
,2), (12 
(20,1) 
(28,2) 
(36,3) 
(44,2) 
(52,1) 
(60,3) 
(68,3) 
(76,3) 
(84,3) 
(92,3) 
(100,3) , 



1) , (4,0) , (5,1) , (6,2) , 
,1), (13,2), (14,3), (15,2) 



(21,2) 
(29,2) 
(37,3) 
(45,2) 
(53,0) 
(61,3) 
(69,4) 
(77,2) 
(85,3) 
(93,4) 
(101,3) 



(22,2) 
(30,3) 
(38,3) 
(46,2) 
(54,1) 
(62,2) 
(70,4) 
(78,3) 
(86,1) 
(94,4) 
, (102 



(23,1) 
(31,4) 
(39,1) 
(47,3) 
(55,2) 
(63,2) 
(71,2) 
(79,4) 
(87,3) 
(95,4) 
4) , (103,4) , 
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(104 


,4) , 


(105 


,4) , 


(106 


,2) , 


(107 


,2) , 


(108 


,2) , 


(109 


,3) , 


(110 


,2 




(111 


,4) , 


(112 


,4) , 


(113 


,4) , 


(114 


,3) , 


(115 


,4) , 


(116 


,2) , 


(117 


,3 




(118 


,4) , 


(119 


,3) , 


(120 


,4) , 


(121 


,3) , 


(122 


,4) , 


(123 


,4) , 


(124 


,4 




(125 


,3) , 


(126 


,3) , 


(127 


,4) , 


(128 


,4) , 


(129 


,3) , 


(130 


,3) , 


(131 


,2 


40 


(132 


,2) , 


(133 


,3) , 


(134 


,3) , 


(135 


,3) , 


(136 


,4) , 


(137 


,4) , 


(138 


,3 




(139 


,4) , 


(140 


,3) , 


(141 


,3) , 


(142 


,3) , 


(143 


,3) : 


(144 


,3) , 


(145 


,3 




(146 


,2) , 


(147 


,4) , 


(148 


,4) , 


(149 


,3) , 


(150 


,3) , 


(151 


,3) , 


(152 


,4 




(153 


,3) , 


(154 


,3) , 


(155 


,2) , 


(156 


,3) , 


(157 


,3) , 


(158 


,3) , 


(159 


,3 




(160 


,3) , 


(161 


,2) , 


(162 


,3) , 


(163 


,3) , 


(164 


,4) , 


(165 


,2) , 


(166 


,3 


45 


(167 


,4) , 


(168 


,3) , 


(169 


,3) , 


(170 


,3) , 


(171 


,3) , 


(172 


,3) , 


(173 


,3 




(174 


,3) , 


(175 


,2) , 


(176 


,1) , 


(177 


,3) , 


(178 


,3) , 


(179 


,3) , 


(180 


,3 




(181 


,3) , 


(182 


,3) , 


(183 


,3) , 


(184 


,3) , 


(185 


,3) , 


(186 


,3) , 


(187 


,3 




(188 


,3) , 


(189 


,2) , 


(190 


,2) , 


(191 


,2) , 


(192 


,3) , 


(193 


,3) , 


(194 


,2 




(195 


,2) , 


(196 


,3) , 


(197 


,2) , 


(198 


,3) , 


(199 


,2) , 


(200 


,2) , 


(201 


,2 


50 


(202 


,3) , 


(203 


,2) , 


(204 


,2) , 


(205 


,2) , 


(206 


,2) , 


(207 


,2) , 


(208 


,1 




(209 


,0) , 


(210 


,1) , 


(211 


,2) , 


(212 


,2) , 


(213 


,1) , 


(214 


,2) , 


(215 


,3 




(216 


,3) , 


(217 


,3) , 


(218 


,2) , 


(219 


,2) , 


(220 


,2) , 


(221 


,3) , 


(222 


,4 




(223 


,2) , 


(224 


,3) , 


(225 


,4) , 


(226 


,4) , 


(227 


,4) , 


(228 


,2) , 


(229 


,3 




(230 


,3) , 


(231 


,3) , 


(232 


,3) , 


(233 


,3) , 


(234 


,2) , 


(235 


,3) , 


(236 


,4 


55 


(237 


,4) , 


(238 


,4) , 


(239 


,3) , 


(240 


,5) , 


(241 


,3) , 


(242 


,3) , 


(243 


,3 




(244 


,4) , 


(245 


,4) , 


(246 


,5) , 


(247 


,5) , 


(248 


,3) , 


(249 


,1) , 


(250 


,3 




(251 


,3) , 


(252 


,2) , 


(253 


,3) , 


(254 


,2) , 


(255 


,3) , 


(256 


,4) , 


(257 


,4 




(258 


,4) , 


(259 


,4) , 


(260 


,4) , 


(261 


,4) , 


(262 


,3) , 


(263 


,4) , 


(264 


,3 




(265 


,4) , 


(266 


,4) , 


(267 


,4) , 


(268 


,4) , 


(269 


,2) , 


(270 


,2) , 


(271 


,3 


60 


(272 


,4) , 


(273 


,4) , 


(274 


,4) , 


(275 


,4) , 


(276 


,5) , 


(277 


,4) , 


(278 


,4 




(279 


,5) , 


(280 


,5) , 


(281 


,4) , 


(282 


,4) , 


(283 


,4) , 


(284 


,5) , 


(285 


,3 




(286 


.6) , 


(287 


,4) , 


(288 


,4) , 


(289 


,4) , 


(290 


,2) , 


(291 


,2) , 


(292 


,3 




(293 


,2) , 


(294 


,4) , 


(295 


,4) , 


(296 


,4) , 


(297 


,3) , 


(298 


,4) , 


(299 


,2 




(300 


,3) , 


(301 


,4) , 


(302 


,3) , 


(303 


,4) , 


(304 


,3) , 


(305 


,4) , 


(306 


,4 


65 


(307 


,4) , 


(308 


,3) , 


(309 


,3) , 


(310 


,4) , 


(311 


,3) , 


(312 


,4) , 


(313 


,3 




(314 


,3) , 


(315 


,4) , 


(316 


,4) , 


(317 


,5) , 


(318 


,4) , 


(319 


,4) , 


(320 


,5 




(321 


,5) , 


(322 


,5) , 


(323 


,4) , 


(324 


,4) , 


(325 


,4) , 


(326 


,3) , 


(327 


,5 




(328 


,5) , 


(329 


,4) , 


(330 


,4) , 


(331 


,5) , 


(332 


,5) , 


(333 


,5) , 


(334 


,5 




(335 


,3) , 


(336 


,5) , 


(337 


,5) , 


(338 


,5) , 


(339 


,4) , 


(340 


,5) , 


(341 


,5 


70 


(342 


,5) , 


(343 


,4) , 


(344 


,4) , 


(345 


,4) , 


(346 


,4) , 


(347 


,3) , 


(348 


,5 




(349 


,5) , 


(350 


,5) , 


(351 


,5) , 


(352 


,3) , 


(353 


,5) , 


(354 


,3) , 


(355 


,2 




(356 


,2) , 


(357 


,3) , 


(358 


,3) , 


(359 


,3) , 


(360 


,4) , 


(361 


,4) , 


(362 


,3 




(363 


,4) , 


(364 


,3) , 


(365 


,3) , 


(366 


,3) , 


(367 


,3) , 


(368 


,3) , 


(369 


,3 




(370 


,2) , 


(371 


,4) , 


(372 


,4) , 


(373 


,3) , 


(374 


,3) , 


(375 


,3) , 


(376 


,3 


75 


(377 


4) , 


(378 


,3) , 


(379 


,4) , 


(380 


,4) , 


(381 


,5) , 


(382 


,3) , 


(383 


,5 




(384 


,4) , 


(385 


,5) , 


(386 


,4) , 


(387 


,3) , 


(388 


,4) , 


(389 


,4) , 


(390 


,5 
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(391 


,5) , 


(392 


,4) , 


(393 


,5) , 


(394 


,5) , 


(395 


,4) , 


(396 


,4) , 


(397 


,3) , 




(398 


,5) , 


(399 


,5) , 


(400 


,5) , 


(401 


,5) , 


(402 


,5) , 


(403 


,4) , 


(404 


,4) , 




(405 


,4) , 


(406 


,4) , 


(407 


,4) , 


(408 


,4) , 


(409 


,4) , 


(410 


,4) , 


(411 


,4) , 


80 


(412 


,5) , 


(413 


,5) , 


(414 


,5) , 


(415 


,4) , 


(416 


,4) , 


(417 


,3) , 


(418 


,3) , 




(419 


,4) , 


(420 


,4) , 


(421 


,5) , 


(422 


,5) , 


(423 


,5) , 


(424 


,4) , 


(425 


,4) , 




(426 


4) , 


(427 


,5) , 


(428 


,4) , 


(429 


,4) , 


(430 


,4) , 


(431 


,4) , 


(432 


,5) , 




(433 


,5) , 


(434 


,5) , 


(435 


,5) , 


(436 


,4) , 


(437 


,5) , 


(438 


,5) , 


(439 


,5) , 




(440 


,4) , 


(441 


,4) , 


(442 


,4) , 


(443 


,4) , 


(444 


,4) , 


(445 


,5) , 


(446 


,5) , 


85 


(447 


,4) , 


(448 


,4) , 


(449 


,4) , 


(450 


,4) , 


(451 


,3) , 


(452 


,2) , 


(453 


,3) , 




(454 


,3) , 


(455 


,3) , 


(456 


,3) , 


(457 


,3) , 


(458 


,2) , 


(459 


,3) , 


(460 


,3) , 




(461 


,4) , 


(462 


,2) , 


(463 


,3) , 


(464 


,4) , 


(465 


,3) , 


(466 


,3) , 


(467 


,3) , 




(468 


,3) , 


(469 


,3) , 


(470 


,3) , 


(471 


,3) , 


(472 


,2) , 


(473 


,3) , 


(474 


,4) , 




(475 


,4) , 


(476 


,4) , 


(477 


,4) , 


(478 


,5) , 


(479 


,5) , 


(480 


,4) , 


(481 


,4) , 


90 


(482 


,5) , 


(483 


,4) , 


(484 


,5) , 


(485 


,4) , 


(486 


,4) , 


(487 


,4) , 


(488 


,5) , 




(489 


,5) , 


(490 


,4) , 


(491 


,4) , 


(492 


,4) , 


(493 


,4) , 


(494 


,4) , 


(495 


,4) , 




(496 


,3) , 


(497 


,5) , 


(498 


,4) , 


(499 


,4) , 


(500 


,4) , 


(501 


,3) , 


(502 


,3) , 




(503 


,4) , 


(504 


,4) , 


(505 


,4) , 


(506 


,4) , 


(507 


,3) , 


(508 


,5) , 


(509 


,5) , 




(510 


,4) , 


(511 


,4) , 


(512 


,5) , 


(513 


,4) , 


(514 


,4) , 


(515 


,5) , 


(516 


,5) , 


95 


(517 


,5) , 


(518 


,5) , 


(519 


,4) , 


(520 


,4) , 


(521 


A) , 


(522 


,4) , 


(523 


,4) , 




(524 


,4) , 


(525 


,3) , 


(526 


,5) , 


(527 


,5) , 


(528 


,5) , 


(529 


,5) , 


(530 


,5) , 




(531 


,4) , 


(532 


,4) , 


(533 


,5) , 


(534 


,4) , 


(535 


,4) , 


(536 


,4) , 


(537 


,4) , 




(538 


,4) , 


(539 


,4) , 


(540 


,4) , 


(541 


,4) , 


(542 


,4) , 


(543 


,4) , 


(544 


,4) , 




(545 


,4) , 


(546 


,4) , 


(547 


,4) , 


(548 


,5) , 


(549 


,4) , 


(550 


,4) , 


(551 


,3) , 


100 


(552 


,4) , 


(553 


,4) , 


(554 


,3) , 


(555 


,4) , 


(556 


,4) , 


(557 


,3) , 


(558 


,3) , 




(559 


,5) , 


(560 


,4) , 


(561 


,5) , 


(562 


,4) , 


(563 


,5) , 


(564 


,4) , 


(565 


,4) , 




(566 


,5) , 


(567 


,4) , 


(568 


,4) , 


(569 


,4) , 


(570 


,3) , 


(571 


,4) , 


(572 


,4) , 




(573 


,4) , 


(574 


,5) , 


(575 


,5) , 


(576 


,4) , 


(577 


,4) , 


(578 


,4) , 


(579 


,4) , 




(580 


,4) , 


(581 


,4) , 


(582 


,2) , 


(583 


,1) , 


(584 


,3) , 


(585 


,3) , 


(586 


,3) , 


105 


(587 


,3) , 


(588 


,3) , 


(589 


,3) , 


(590 


,3) , 


(591 


,3) , 


(592 


,3) , 


(593 


,3) , 




(594 


,3) , 


(595 


,3) , 


(596 


,2) , 


(597 


,2) , 


(598 


,4) , 


(599 


,4) , 


(600 


,4) , 




(601 


,4) , 


(602 


,4) , 


(603 


,4) , 


(604 


,4) , 


(605 


,4) , 


(606 


,4) , 


(607 


,4) , 




(608 


,4) , 


(609 


,4) , 


(610 


,3) , 


(611 


,3) , 


(612 


,3) , 


(613 


,4) , 


(614 


,4) , 




(615 


,3) , 


(616 


,4) , 


(617 


,4) , 


(618 


,5) , 


(619 


,3) , 


(620 


,4) , 


(621 


,4) , 


110 


(622 


,5) , 


(623 


,5) , 


(624 


,4) , 


(625 


,4) , 


(626 


,4) , 


(627 


,4) , 


(628 


,4) , 




(629 


,4) , 


(630 


,3) , 


(631 


,4) , 


(632 


,5) , 


(633 


,4) , 


(634 


,4) , 


(635 


,4) , 




(636 


,4) , 


(637 


,4) , 


(638 


,4) , 


(639 


,4) , 


(640 


,5) , 


(641 


,4) , 


(642 


,4) , 




(643 


,4) , 


(644 


,3) , 


(645 


,4) , 


(646 


,5) , 


(647 


,4) , 


(648 


,4) , 


(649 


,4) , 




(650 


,4) , 


(651 


,4) , 


(652 


,4) , 


(653 


,4) , 


(654 


,4) , 


(655 


,4) , 


(656 


,4) , 


115 


(657 


,4) , 


(658 


,4) , 


(659 


,4) , 


(660 


,4) , 


(661 


,4) , 


(662 


,4) , 


(663 


,4) , 




(664 


4) , 


(665 


,4) , 


(666 


,4) , 


(667 


,3) , 


(668 


,3) , 


(669 


,3) , 


(670 


,2) , 




(671 


,4) , 


(672 


,4) , 


(673 


,4) , 


(674 


,4) , 


(675 


,4) , 


(676 


,4) , 


(677 


,4) , 
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(678 


,3) , 


(679 


,4) , 


(680 


,4) , 


(681 


,3) , 


(682 


,5) , 


(683 


,5) , 


(684 


,4) , 




(685 


,4) , 


(686 


,3) , 


(687 


,3) , 


(688 


,4) , 


(689 


,3) , 


(690 


,5) , 


(691 


,4) , 


120 


(692 


,4) , 


(693 


,4) , 


(694 


,5) , 


(695 


,4) , 


(696 


,4) , 


(697 


,4) , 


(698 


,4) , 




(699 


,4) , 


(700 


,4) , 


(701 


,4) , 


(702 


,4) , 


(703 


,4) , 


(704 


,4) , 


(705 


,4) , 




(706 


,4) , 


(707 


,4) , 


(708 


,4) , 


(709 


,4) , 


(710 


,5) , 


(711 


,4) , 


(712 


,4) , 




(713 


,4) , 


(714 


,4) , 


(715 


,4) , 


(716 


,4) , 


(717 


,4) , 


(718 


,4) , 


(719 


,4) , 




(720 


,4) , 


(721 


,3) , 


(722 


,4) , 


(723 


,4) , 


(724 


,4) , 


(725 


,3) , 


(726 


,3) , 


125 


(727 


,4) , 


(728 


,4) , 


(729 


,4) , 


(730 


,3) , 


(731 


,2) , 


(732 


,3) , 


(733 


,3) , 




(734 


,2) , 


(735 


,2) , 


(736 


,3) , 


(737 


,2) , 


(738 


,3) , 


(739 


,2) , 


(740 


,4) , 




(741 


,4) , 


(742 


,4) , 


(743 


,3) , 


(744 


,4) , 


(745 


,2) , 


(746 


,4) , 


(747 


,4) , 




(748 


.4) , 


(749 


,4) , 


(750 


,4) , 


(751 


,4) , 


(752 


,4) , 


(753 


,4) , 


(754 


,4) , 




(755 


,4) , 


(756 


,4) , 


(757 


,4) , 


(758 


,4) , 


(759 


,4) , 


(760 


,4) , 


(761 


,4) , 


130 


(762 


,4) , 


(763 


,4) , 


(764 


,4) , 


(765 


,4) , 


(766 


,4) , 


(767 


,4) , 


(768 


,4) , 




(769 


,4) , 


(770 


,4) , 


(771 


,4) , 


(772 


,4) , 


(773 


,4) , 


(774 


,4) , 


(775 


,4) , 




(776 


,4) , 


(777 


,4) , 


(778 


,4) , 


(779 


,4) , 


(780 


,3) , 


(781 


,3) , 


(782 


,4) , 




(783 


,4) , 


(784 


,4) , 


(785 


,4) , 


(786 


,3) , 


(787 


,3) , 


(788 


,4) , 


(789 


,3) , 




(790 


,2) , 


(791 


,3) , 


(792 


,4) , 


(793 


,4) , 


(794 


,4) , 


(795 


,4) , 


(796 


,3) , 


135 


(797 


,4) , 


(798 


,4) , 


(799 


,4) , 


(800 


,4) , 


(801 


,3) , 


(802 


,4) , 


(803 


,4) , 




(804 


,4) , 


(805 


,4) , 


(806 


,4) , 


(807 


,4) , 


(808 


,4) , 


(809 


,3) , 


(810 


A) , 




(811 


,4) , 


(812 


,3) , 


(813 


,5) , 


(814 


,3) , 


(815 


,3) , 


(816 


,4) , 


(817 


,4) , 




(818 


,4) , 


(819 


,4) , 


(820 


,4) , 


(821 


,5) , 


(822 


,4) , 


(823 


,4) , 


(824 


,4) , 




(825 


,4) , 


(826 


,4) , 


(827 


,4) , 


(828 


,4) , 


(829 


,3) , 


(830 


,4) , 


(831 


,3) , 


140 


(832 


.3) , 


(833 


,4) , 


(834 


,4) , 


(835 


,4) , 


(836 


,4) , 


(837 


,4) , 


(838 


,4) , 




(839 


,4) , 


(840 


,5) , 


(841 


,4) , 


(842 


,4) , 


(843 


,4) , 


(844 


,4) , 


(845 


,4) , 




(846 


,3) , 


(847 


,4) , 


(848 


,4) , 


(849 


,4) , 


(850 


,4) , 


(851 


,3) , 


(852 


,4) , 




(853 


,4) , 


(854 


,4) , 


(855 


,4) , 


(856 


,3) , 


(857 


,4) , 


(858 


,4) , 


(859 


,4) , 




(860 


,4) , 


(861 


,4) , 


(862 


.4) , 


(863 


,3) , 


(864 


,4) , 


(865 


,3) , 


(866 


,4) , 


145 


(867 


,4) , 


(868 


,4) , 


(869 


,3) , 


(870 


,4) , 


(871 


,4) , 


(872 


,3) , 


(873 


,3) , 




(874 


,3) , 


(875 


,4) , 


(876 


,3) , 


(877 


,4) , 


(878 


,3) , 


(879 


,2) , 


(880 


,2) , 




(881 


,3) , 


(882 


,2) , 


(883 


,2) , 


(884 


,3) , 


(885 


,3) , 


(886 


,4) , 


(887 


,4) , 




(888 


,4) , 


(889 


,4) , 


(890 


,4) , 


(891 


,4) , 


(892 


,4) , 


(893 


,3) , 


(894 


,3) , 




(895 


,3) , 


(896 


,3) , 


(897 


,4) , 


(898 


,3) , 


(899 


,4) , 


(900 


,4) , 


(901 


,3) , 


150 


(902 


,3) , 


(903 


,4) , 


(904 


,4) , 


(905 


,4) , 


(906 


,3) , 


(907 


,4) , 


(908 


,4) , 




(909 


,3) , 


(910 


,4) , 


(911 


,3) , 


(912 


,3) , 


(913 


,3) , 


(914 


,4) , 


(915 


,4) , 




(916 


,4) , 


(917 


,4) , 


(918 


,3) , 


(919 


,4) , 


(920 


,4) , 


(921 


,4) , 


(922 


,4) , 




(923 


,4) , 


(924 


,3) , 


(925 


,3) , 


(926 


,4) , 


(927 


,4) , 


(928 


,4) , 


(929 


,4) , 




(930 


,4) , 


(931 


,3) , 


(932 


,3) , 


(933 


,4) , 


(934 


,4) , 


(935 


,4) , 


(936 


,4) , 


155 


(937 


,4) , 


(938 


,4) , 


(939 


,4) , 


(940 


,4) , 


(941 


,4) , 


(942 


,4) , 


(943 


,4) , 




(944 


,3) , 


(945 


,3) , 


(946 


,4) , 


(947 


,3) , 


(948 


,3) , 


(949 


,3) , 


(950 


,4) , 




(951 


4) , 


(952 


,4) , 


(953 


,3) , 


(954 


,4) , 


(955 


,4) , 


(956 


,3) , 


(957 


,3) , 




(958 


,3) , 


(959 


,4) , 


(960 


,4) , 


(961 


,4) , 


(962 


,4) , 


(963 


,4) , 


(964 


,4) , 
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(965 ,4) , 


(966 ,4) 




(967,4) 




;968 ,4) , 


(969 ,4) , 


(970 


,3) , (971, 


4) , 




160 


(972,4) , 


(973 ,3) 




(974 ,4) 




;975 ,3) , 


(976 ,4) , 


(977 ,3) , (978 , 


3) , 






(979 ,4) , 


(980 ,4) 




(981 ,4) 




;982 ,4) , 


(983 


,3), 


(984 ,3) , (985 , 


4) , 






(986 ,4) , 


(987,3) 




(988 ,3) 




;989 ,4) , 


(990 ,3) , 


(991 


,3) , (992 , 


4) , 






(993 ,4) , 


(994 ,3) 




(995 ,3) 




;996 ,3) , 


(997 ,4) , 


(998 


,4) , (999 , 


4) , 






(1000 ,3) 


, (1001 


3) 


, (1002 


,3) 


, (1003 


:3) , 


;i004 


,3) , 


(1005 ,3) , 


(1006 


,4 


165 


(1007 ,2) 


, (1008 


4) 


, (1009 


,2) 


, (1010 


,2) , 


;ioii 


,2) , 


(1012 ,3) , 


(1013 


,3 




(1014 ,3) 


, (1015 


4) 


, (1016 


,4) 


, (1017 


,3) , 


;i018 


,3) , 


(1019 ,3) , 


(1020 


,3 




(1021 ,4) 


, (1022 


3) 


, (1023 


,3) 


, (1024 


,3) , 


^1025 


,4) , 


(1026 ,4) , 


(1027 


,4 




(1028 ,3) 


, (1029 


4) 


, (1030 


,4) 


, (1031 


,4) , 


^1032 


,2) , 


(1033 ,4) , 


(1034 


,3 




(1035 ,3) 


, (1036 


3) 


, (1037 


>3) 


, (1038 


,3) , 


;i039 


,4) , 


(1040 ,3) , 


(1041 


,4 


170 


(1042 ,3) 


, (1043 


4) 


, (1044 


>3) 


, (1045 


,3) , 


;i046 


,4) , 


(1047 ,4) , 


(1048 


,4 




(1049 ,3) 


, (1050 


4) 


, (1051 


,4) 


, (1052 


,4) , 


;i053 


,4) , 


(1054 ,4) , 


(1055 


,4 




(1056 ,3) 


, (1057 


3) 


, (1058 


,4) 


, (1059 


,4) , 


^1060 


,3) , 


(1061 ,4) , 


(1062 


,3 




(1063 ,3) 


, (1064 


3) 


, (1065 


,4) 


, (1066 


,3) , 


^1067 


,3) , 


(1068 ,3) , 


(1069 


,4 




(1070 ,3) 


, (1071 


4) 


, (1072 


>3) 


, (1073 


,3) , 


^1074 


,3) , 


(1075 ,3) , 


(1076 


,3 


175 


(1077 ,3) 


, (1078 


4) 


, (1079 


,3) 


, (1080 


,4) , 


^1081 


,3) , 


(1082 ,4) , 


(1083 


,4 




(1084 ,3) 


, (1085 


3) 


, (1086 


,3) 


, (1087 


,3) , 


^1088 


,2) , 


(1089 ,4) , 


(1090 


,3 




(1091 ,4) 


, (1092 


3) 


, (1093 


,4) 


, (1094 


,3) , 


^1095 


,3) , 


(1096 ,3) , 


(1097 


,4 




(1098 ,3) 


, (1099 


3) 


, (1100 


>3) 


, (1101 


,4) , 


^1102 


,3) , 


(1103 ,3) , 


(1104 


,3 




(1105 ,3) 


, (1106 


3) 


, (1107 


,2) 


, (1108 


,3) , 


^1109 


,3) , 


(1110 ,3) , 


(1111 


,3 


180 


(1112 ,3) 


, (1113 


3) 


, (1114 


,3) 


, (1115 


,3) , 


^1116 


,3) , 


(1117 ,4) , 


(1118 


,4 




(1119 ,3) 


, (1120 


3) 


, (1121 


,4) 


, (1122 


,3) , 


^1123 


,3) , 


(1124 ,3) , 


(1125 


,3 




(1126 ,3) 


, (1127 


4) 


, (1128 


>3) 


, (1129 


,3) , 


^1130 


,3) , 


(1131 ,3) , 


(1132 


,3 




(1133 ,3) 


, (1134 


3) 


, (1135 


,3) 


, (1136 


,3) , 


^1137 


,3) , 


(1138 ,3) , 


(1139 


,3 




(1140 ,2) 


, (1141, 


4) 


, (1142 


,4) 


, (1143 


,3) , 


^1144 


,3) , 


(1145 ,4) , 


(1146 


,3 


185 


(1147 ,3) 


, (1148 


3) 


, (1149 


.3) 


, (1150 


,4) , 


^1151 


,3) , 


(1152 ,3) , 


(1153 


,3 




(1154 ,4) 


, (1155 


3) 


, (1156 


,3) 


, (1157 


,3) , 


;il58 


,3) , 


(1159 ,3) , 


(1160 


,4 




(1161 ,3) 


, (1162 


3) 


, (1163 


,3) 


, (1164 


,2) , 


^1165 


,3) , 


(1166 ,3) , 


(1167 


,3 




(1168 ,3) 


, (1169 


3) 


, (1170 


,3) 


, (1171 


,2) , 


^1172 


,1) , 


(1173 ,3) , 


(1174 


,3 




(1175 ,3) 


, (1176 


3) 


, (1177 


>3) 


, (1178 


,3) , 


^1179 


,3) , 


(1180 ,3) , 


(1181 


,3 


190 


(1182 ,3) 


, (1183 


3) 


, (1184 


,2) 


, (1185 


,3) , 


;il86 


,3) , 


(1187 ,4) , 


(1188 


,2 




(1189 ,3) 


, (1190 


4) 


, (1191 


,3) 


, (1192 


,3) , 


^1193 


,3) , 


(1194 ,3) , 


(1195 


,3 




(1196 ,3) 


, (1197 


3) 


, (1198 


,3) 


, (1199 


,3) , 


^1200 


,3) , 


(1201 ,3) , 


(1202 


,3 




(1203 ,3) 


, (1204 


3) 


, (1205 


>3) 


, (1206 


,2) , 


^1207 


,3) , 


(1208 ,2) , 


(1209 


,3 




(1210 ,3) 


, (1211 


2) 


, (1212 


>3) 


, (1213 


,2) , 


^1214 


,3) , 


(1215 ,3) , 


(1216 


,2 


195 


(1217 ,3) 


, (1218 


3) 


, (1219 


,3) 


, (1220 


,3) , 


^1221 


,3) , 


(1222 ,3) , 


(1223 


,3 




(1224 ,3) 


, (1225 


3) 


, (1226 


,2) 


, (1227 


,3) , 


^1228 


,3) , 


(1229 ,2) , 


(1230 


,2 




(1231 ,3) 


, (1232 


3) 


, (1233 


,2) 


, (1234 


,2) , 


^1235 


,3) , 


(1236 ,3) , 


(1237 


,2 




(1238 ,2) 


, (1239 


3) 


, (1240 


,2) 


, (1241 


,3) , 


;i242 


,2) , 


(1243 ,2) , 


(1244 


,2 




(1245 ,2) 


, (1246 


3) 


, (1247 


,2) 


, (1248 


,2) , 


^1249 


,2) , 


(1250 ,2) , 


(1251 


,2 
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200 (1252,1)] 



///////////////////////////////////////////////////////////////^^^ 

# The atlas of graphs is ordered by number of vertices , then number 

^of 

205 # edges. To make our search for minimum ranks more efficient , we 

# store the first and last indices of the graphs having a certain 

# number of vertices and edges. 

# graph_help [n ] [m] is a tuple (i,j), meaning that the graphs with n 
210 # vertices and m edges start at index i and end on index j in the 

# atlas of graphs. 

///////////////////////////m^^^^ 



215 



graph_help = { 
1: {0: (1,1)}, 

2: {0: (2,2) , 1: (3,3)}, 

3: {0: (4,4) , 1: (5,5) , 2: (6,6) , 3: (7,7)}, 

4: {0: (8,8), 1: (9,9), 2: (10,11), 3: (12,14), 4: (15,16), 

5: (17,17) , 6: (18,18)}, 

5: {0: (19,19), 1: (20,20), 2: (21,22), 3: (23,26), 4: (27,32), 5: 

(33,38), 6: (39,44), 7: (45,48), 8: (49,50), 9: (51,51), 10: 
(52,52)}, 

6: {0: (53,53), 1:(54,54), 2:(55,56), 3: (57,61), 4: (62,70), 
5: (71,85), 6: (86,106), 7: (107,130), 8: (131,154), 9: ^ 
-^(155 ,175) , 

10: (176,190), 11: (191,199), 12: (200,204), 13: (205,206), 
14: (207 ,207) , 15: (208 ,208) } , 

7: {0: (209,209), 1: (210,210), 2: (211,212), 3: (213,217), 

4: (218,227), 5: (228,248), 6: (249,289), 7: (290,354), 8: ^ 

■-^(355 ,451) , 

9: (452,582), 10: (583,730), 11: (731,878), 12: (879,1009), 
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13: (1010,1106), 14: (1107,1171), 15: (1172,1212), 16: 

-^(1213 ,1233) , 
17: (1234,1243), 18: (1244,1248), 19: (1249,1250), 20: 

■-^(1251,1251) , 
21: (1252,1252)}} 



def get_mr_from_list (graph) : 

?5 !J 

245 Return the minimum rank of a graph of order 7 or less from the 

list min_ranks . 



INPUT: 

graph — the graph whose minimum rank is to be 
260 found 



OUTPUT: 

the minimum rank if the graph is in the list 
or False if it is not 

EXAMPLES: 

sage: get_mr_from_list (Graph ( { : [2 ,3],1:[2],2:[3,4],3:[4]})) 
3 

sage : get_mr_from_list( graphs. PathGraph (8) ) 
False 



•)■) !? 5? 



#check to make sure graph can he found in list 
265 if graph . order >7: 

return False 



order = graph . order ( ) ^umher of vertices 
size = graph, size #number of edges 

270 

starting_index , ending_index = graph_help [ order ][ size ] 
import pdb 

#look for graph and return minimum rank 
for i in [ starting_index . . ending_indcx ] : 
276 if graph . is_isomorphic ( atlas -graphs [ i ]) : 

return min_ranks [ i ] [ 1 ] 
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raise ValueError , "This should never happen!" 



280 

def zerosgame ( graph , i n i t i a 1 _s e t = [] ) : 

5? 5? 

Apply the color— change rule to a given graph given an optional 
initial set . 

285 

INPUT: 

graph — the graph on which to apply the rule 
initial_sct — the set of "zero" (black) vertices in the <- 
■^graph 

290 OUTPUT: 

the list of zero (black) vertices in the resulting derived 
coloring 

EXAMPLES: 

295 sage: zerosgame ( graphs . PathGraph (5) ) 

[] 

sage: zerosgame ( graphs . PathGraph (5) ,[0]) 
[0, 1, 2, 3, 4] 

?? 55 55 

300 

new_zero_set=set ( initial_set ) 
zero_set=set ( [] ) 
zer o_neighbor s={} 
active_zero_set = set([]) 
305 inactive_zero_set = set([]) 

another_r un=True 
while another_run : 

another_run=False 

# Add the new zero vertices 
310 zero_set . update ( new_zero_set ) 

active_zero_set.update(new_zero_set) 

active_zero_set.difference_update(inactive_zero_set) 
zero_neighbor s . update ( [ [ i , 

set (graph, neighbors ( i ) ) . difference (zero_set ) ] 
315 for i in new_zero_set ] ) 

# Find the next set of zero vertices 
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new_zero_set . clear () 
inactive_zero_set . clear () 
for V in active_zero_set : 
320 zero_neighbors[v]. difference .update (zero_set) 

if len ( zero.neighbors [v] ) ==1: 

new_zero_set . add (zero.neighbors [v] . pop () ) 
inactive_zero_set .add(v) 
another_run=True 
326 return 1 i s t ( zer o _set ) 

def find_zero_forcing_set (graph , bound=None) : 

5? ?! !5 

330 Return a zero forcing set of minimum order that also has order 
less than the given bound. 

INPUT: 

graph — the graph on which to find the zero — forcing set 
336 bound — the maximum acceptable order for a zero — forcing set 

OUTPUT: 

a zero— forcing set of minimum order that also has order less 
than the bound if one exists ; False if no such zero — forcing set 
340 can be found 

EXAMPLES: 

sage : find_zero_forcing_set (graphs. CompleteGraph (5) ) 
{0, 1, 2, 3} 

346 sage: find _zero _f or cing_set ( graphs . CompleteGraph ( 5 ), 2 ) 

False 

7? V V 

order=graph . order () 
if bound is None: 
360 bound = order 

if bound<0: 
bound=l 
vertices=graph . vertices () 
mindegree=min ( graph . degree ( ) ) 
366 for i in [ mindegrec .. bound ] : 

for subset in Subsets ( vertices , i ) : 
outcome=zerosgame ( graph , subset ) 
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360 



if len ( outcome )=order : 
return subset 

return False 



def find_Z ( graph) : 

17 -7 71 

365 Returns the order of a smallest zero— forcing set of a graph 

INPUT: 

graph — the graph on which to find a smallest zero — forcinj 

■^set 

370 OUTPUT: 

the minimum possible order for a zero— forcing set of the <- 
'— >graph 

EXAMPLES: 

sage : find_Z (graphs . CompleteGraph (5) ) 

376 4 
?! 55 5? 

return len (find_zero_forcing_set( graph)) 



380 



def has_forbidden_induced_subgraph ( graph) : 

5? 5? 55 

Check for a forbidden induced subgraph (a path on 4 vertices , 
fish , dart , or complete tripartite graph) . 

INPUT: 

graph — the graph to be checked 
OUTPUT: 

True if the graph contains an induced copy of P-4 , fish , dart 

^ ^ 5 

or K_{3,3,3}; False if it does not 
EXAAIPLES: 

sage : has_for bidden_induced_subgr aph (graphs . CompleteGraph (10) 

12 



395 False 

sage : has_forbidden_induced_subgr aph (graphs . PathGraph (3) ) 

False 

sage: K333 = Graph ({ : [3,4,5,6,7,8], 1: [3,4,5,6,7,8], 2: 
-^[3,4,5,6,7,8], 3: [6,7,8], 4: [6,7,8], 5: [6,7,8]}) 





sage : 


has_forbidden. 


.induced_subgraph (K333) 


400 


True 








sage : 


g = Graph ( { : 


[1,2,3,4], 1: [2]}) # fi 




sage : 


has_forbidden. 


.induced_subgraph(g) 




True 








sage : 


g . add_cdge ( (0 


,6)) 


405 


sage : 


has_for bidden. 


.induced_subgraph(g) 




True 







?? •>■) 5? 



order = graph . or der ( ) 
vertices = graph . vertices () 
410 path = atlas_graphs [14] 

fish = atlas_graphs [34] 
dart = atlas_graphs [40] 

K333 = Graph ({0: [3,4,5,6,7,8], 1: [3,4,5,6,7,8], \ 
2: [3,4,5,6,7,8] , 3: [6,7,8] , \ 
4: [6,7,8] , 5: [6,7,8]}) 

if order < 4: 

return False 
for sub_vertices in Combinations ( vertices ,4) : 

# Finds all order 4 induced subgraphs 

420 if graph . subgraph ( sub_vert ices ). is_isomorphic ( path ) : 

return True 
if order < 5: 

return False 
for sub_vcrticcs in Combinations ( vertices ,5) : 
425 # Finds all order 5 induced subgraphs 

if graph. subgraph(sub .vertices ) . is_isomorphic (dart) : 
return True 

if gr aph. subgraph(sub .vertices) .is_isomorphic(fish) : 
return True 
430 if order < 9: 

return False 

for sub_vcrticcs in Combinations ( vertices ,9) : 

# Finds all order 9 induced subgraphs 

if graph, subgraph (sub_vertices ) . is_isomorphic (K333) : 
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435 return True 

return False 



def is_outerplanar (graph) : 

;7 ;7 1? 

Check if the graph is outer— planar 
INPUT: 

graph — is the graph to be checked 
OUTPUT: 

True if the graph is outer —planar ; False if it is not 
EXAAIPLES: 

sage: is.out er planar (Graph ({ :[ 1 ,2,3,4],1:[2,4],2:[3],3:[4]}) 
False 

sage : is_outerplanar (graphs . CompleteGraph (3) ) 
True 

?! !5 ;? 

h = graph . copy () 

h. delete_vertices ( [v for v in h.vertices() if h.degree(v) = 0]) 
if h . order ( ) ==0: 

return True 
h.set_boundary(h. vertices ()) 
return h. is_circular _planar ( ordered=False ) 



def find_edge_max_clique (edge , graph) : 

55 5! ?) 

Given an edge and a graph , return a maximal clique of the graph 
that contains the edge. 

INPUT: 

edge — the edge for which to find the maximal clique 
graph — the graph containing the edge 

OUTPUT: 

the list of vertices in a maximal clique that contains the < 
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475 



^edge 

if the edge is not in the graph, it will return None 



EXAMPLES: 

sage: find_edge_max_clique ( ( 1 , 2 ) , graphs . CompleteGraph ( 5 ) ) 
[0, 1, 2, 3, 4] 

sage: find_edge_max_clique ( ( 1 ,3) , graphs . PathGraph (5) ) 

?5 5! 

vertexl=edge [0] # first vertex of edge 

vertex2=edge [ 1 ] # second vertex of edge 

pot _cliques=graph . cliques .containing .vertex (vert ex 1 ) 

# sort the cliques containing vertexl by order, largest first 
pot_cliques . sort (key=len , rever se=True ) 

for clique in pot_cliques: 
if vertex2 in clique : 
return clique 
return None 



495 def edge_clique_cover (graph , bound=None) : 

?5 ?5 71 

Assuming the graph is connected , return an edge clique cover for 
the graph if the number of covering cliques is at most bound; 
otherwise, returns None. 

600 

INPUT: 

graph — the graph 

bound — the maximum number of cliques to consider 
605 OUTPUT: 

If a clique cover is found that has at most bound cliques , •s-^ 
-^•the 

clique cover is returned as a list of lists , each sublist 
being the vertices of a clique . 

610 If a clique cover from this function requires more than bound 

cliques , None is returned. 

EXAMPLES: 
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sage : edge_clique_cover( graphs. PathGraph (3) ) 
[[0, 1], [1, 2]] 

sage : edgc_clique_cover (graphs . CompleteGraph (5) ) 
[[0, 1, 2, 3, 4]] 

sage : edge_clique_cover( graphs. HouseGraph ( ) ) 
[[2, 3, 4], [0, 1], [0, 2], [1, 3]] 

sage: edge_clique_cover (graphs . PetersenGraph ,bound=4) 
sage : 



?5 !J 



# Take care of trivial case 
if graph, size = 0: 
return [] 

max_cliques=graph . cliques () 
max_cliques . sort (key=len) 

largest_clique_vertices = len ( max_cliques [ — 1] ) 
max_cliques = [ sorted ( clique ) for clique in max_cliques 
largest _clique_edges = largest_clique_vertices \ 

*( largest .clique _v ert ices — 1)/2 
edges_of_gr aph=graph . edges ( labels=False ) 
num_edges = graph, size () 

mandatory_cliques = [] 



for V in graph . vertices () : 

# If V is contained in only one clique , then that clique must 

# be in the clique cover 

cliques_containing_v = [c for c in max_cliques if v in c] 
if len ( cliques_containing_v )==1 \ 

and ( cliques_containing_v [0] not in mandatory_cliques<— > 
^): 

mand at or y .clique s . append ( cliques_containing_v [0]) 
for e in graph . edges () : 

# If e is contained in only one clique , then that clique must 

# be in the clique cover 

cliques_containing_e = [c for c in max_cliques 

if e[0] in c and e[l] in c] 

if len ( cliques_containing_e )==1 \ 

and ( cliques_containing_e [0] not in mandatory.cliques*-^ 
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mandatory .cliques . append ( cliques_containing_e [0]) 

# Check to see if mandatory ^cliques contains a clique cover 
edges_in_set_of_cliques = set([]) 
for clique in mandatory_cIiques : 

edges_in_clique = [(clique[i], clique[j]) 

for i in xrange ( len ( clique ) ) 
for j in xrange ( i +1, len ( clique )) ] 
edges_in_set_of_cliques .update(set (edges_in_clique)) 
if len ( edges_in_set _of .cliques ) = num_edges : 

if bound is None or len ( mandatory_cliques ) <= bound: 

return mandatory_cliques 
else : 

# There are too many cliques. Return None to be 

# consistent with the documentation , even though we 

# actually know the clique cover number (and it is <^ 

^greater 

# than bound ) . 
return None 

max_cliques = [c for c in max_cliques if c not in 

^mandatory .cliques ] 
i f bound=None : 

stopping_point=len( max .cliques) 
else : 

stopping_point=min ( len ( max .cliques ) , bound— len 
■^mandatory.cliques ) ) 

starting_point = max(l , ceil (num_edges / largest _clique_edges ) \ 

— len ( mandatory_cliques ) ) 
for i in [ starting_point . . stopping_point ] : 

for set _of .cliques in Combinations ( max_cliques , i ) : 
edges_in_set_of_cliques = set([]) 
for clique in s e t _o f _c 1 i qu e s+mandatory_cliques : 
edges_in_clique = [(clique[i], clique[j]) 

for i in xrange ( len ( clique ) ) 
for j in xrange ( i +1 , len ( clique )) ] 
cdgcs_in_sct_of_cliques.updatc(sct(cdgcs_in_clique)) 
if len ( e d ge s _i n _se t _of _c 1 i q u e s ) = num_edges : 
return set _of _cliques+mandatory_cliques 
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return None 



def edge_clique_cover_approximate (graph , bound=None) : 

55 55 5? 

695 

Returns a decent (though not necessarily tight) upper bound for 
the clique cover number, which is the minimum number of cliques 
necessary to cover all of the edges in a graph 

600 INPUT : 

graph — the graph to be examined 
bound — this argument is ignored 

OUTPUT: 

605 a list of lists of vertices for each clique in the cover 



EXAMPLE: 

sage: sorted ([ sorted ( e ) for e in 

■^edge.clique.cover.approximate (graphs . PathGraph (4) ) ] ) 
[[0, 1], [1, 2], [2, 3]] 

55 55 5! 

610 

free .edges =graph . edges ( labels=False ) 

vertices=graph . vertices () 

vertices . sort (key=graph . degree ) # sort the vertices by degree 
clique_vertices = [] #list of clique vertices 
615 while len ( free_edges ) >0: 

clique_edge=None 
while clique_edge is None: 

# find an edge adjacent to a vertex of minimum degree 
v=vertices [0] 

620 edge_exists , edge=exists ( free_edges , lambda edge: v in •e-^ 

^edge ) 
if edge_exists : 

clique _edge=edge 
else : 

vertices . pop (0) 

625 

edge_clique=find_edge_max_clique(clique_edge , graph) 

clique_vcrtices.append(edgc_clique) 

edges_in_clique = graph . subgraph ( edge_clique ). edges ( labels=<-^ 
-^False) 
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630 



free_cdgcs= [e for e in free_edges if e not in <^ 
^edges_in_clique] 
return clique_vertices 



def min_rank_by_bounds ( graph , tests = [ ' precomputed ' , 'order', 'zero 
^forcing', 'not path', 'forbidden minrank 2', 'not planar', 'not 
router planar', 'clique cover', 'diameter']): 

?5 !J 5? 

Return dictionaries giving the upper and lower bounds from -e-^ 
^running 

636 the specified tests . If tests is not set , then all applicable 

tests are run. 



INPUT: 

graph — the graph for which to find bounds 

640 

OUTPUT: 

a list of 2 dictionaries; the upper and lower bounds, 
^respectively . 



EXAMPLE: 

646 sage: g = Graph({0: [1,2,4,6,7], 1: [3,5,6,7,8], 2: [4,6,8], 

-^3: [4,7,6] , 4: [6] , 5: [6,8,7]}) 
sage: min_rank_by_bounds (g) 
({'zero forcing ': 4}, 
{ ' clique cover ' : 9 , 
'not outer planar ': 6, 
660 ' not path ' : 7 , 

' not planar ' : 5 , 
'order ': 8}) 

sage: min_rank_by_bounds (g , tests=['zero forcing', 'order', 

■^not path ' ] ) 
({'zero forcing': 4}, {'not path': 7, 'order': 8}) 

55 5! ?? 

666 

if isinstance ( tests , str): 
tests = [ tests ] 



order = graph . order ( ) 

lower_bound = {} 
upper_bound = {} 
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' precomputed ' in tests: 
mr = get_mr_from_list (graph) 
i f mr is not False : 

lower_bound [' precomputed ' ] = mr 

upper_bound [ 'precomputed'] = mr 

' order ' in tests : 

upper_bound [' order ' ] = order — 1 
'zero forcing' in tests: 

lower_bound [ ' zero forcing'] = order — find_Z ( graph ) 

# Check if graph is a tree . 

# If yes, then the ZFS will determine minimum rank. 
if graph. is_tree () : 

upper_bound [ ' zero forcing (tree)'] = lower_bound [ ' zero 
^forcing ' ] 

'not path' in tests: 
if graph . diameter ( ) < order — 1: 

upper_bound [ ' not path'] = order — 2 

'forbidden minrank 2' in tests: 
if has_for bidden_induced_subgr aph ( graph ) : 

lower_bound [' forbidden minrank 2'] =3 
else : 

upper_bound [' forbidden minrank 2'] =2 
'diameter' in tests: 

lower_bound [' diameter ' ] = graph . diameter ( ) 
'not planar' in tests: 

# Old versions of Sage assume that planar testing does not 

# have vertices of degree zero. We can delete vertices of 

# degree zero without affecting the planarity . 
h = graph . copy ( ) 

h. delete_vertices ( [v for v in h.vertices() if h.degree(v) = 
--0]) 

if h.order()>0 and h . is _planar ( ) is False: 
upper_bound [ ' not planar'] = order — 4 
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if 'not outer planar' in tests: 

if is_outerplanar (graph) is False: 

upper_bound [ ' not outer planar'] = order — 3 

705 

if 'clique cover' in tests: 

upper_bound [' clique cover'] = len ( edge_clique_cover ( graph ) ) 

return ( lower _bound , upper_bound) 



def find_cut_vertex (graph) : 

)5 5? 5? 

Return a "good" cut— vertex for a graph if it exists; otherwise, 
returns False . 

INPUT: 

graph — the graph on which to find a cut— vertex 
OUTPUT: 

a cut— vertex (if one exists) that either results in 
^components 

of order less than 7 or a minimum of the maximum component 
order ; otherwise False 

EXAMPLES: 

sage : find_cut .vertex (graphs. PathGraph (3) ) 
1 

sage : find_cut_vertex( graphs. PathGraph (20) ) 
9 

sage: [ find_cut_vertex ( graphs . PathGraph ( i ) ) for i in [1..20]] 
[False, False, 1, 1, 1, 1, 1, 1, 2, 3, 4, 5, 6, 6, 7, 7, 8, ^ 
-^8, 9, 9] 

sage : find_cut_vertex (graphs . CompleteGraph (3) ) 
False 



vertices=graph . vertices () 

graph_cc_num=graph . connected_components_number () 
gr aph_order=graph . order () 

#this will hold the "best" cut— vertex and the order of the 
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"-^ lar g est 

^connected component after deletion 
best_v=(False , graph_order ) 

#checks each vertex and determines the best one 
746 for V in vertices: 

g=graph . copy 
g. delete_vertex(v) 
g_cc = g . connected-components 
if len ( g_cc )>graph_cc_num : 
750 # We have a cut—vertex 

max_order = max(len(c) for c in g_cc) 



if max_order <7: 
return v 

766 if max_ordcr<best_v [ 1 ] : 

best_v = (v , max_order ) 



return best_v [0] 



def find_rank_spread ( vertex , graph): 

;7 57 1? 

Returns the exact rank spread for a graph and a vertex (i.e., 
mr (G)— mr(G— v) ) if the minimum ranks of both the graph and the 
766 graph without the vertex can be calculated using the ^ 

■^minrank.bounds 
function . 



INPUT: 

vertex — the vertex 
graph — the graph 



OUTPUT: 

the rank spread and mr( graph— vertex ) if both can be 
^calculated 

using the minrank_bounds program, or False and False if 
■^either cannot 
776 be calculated exactly . 



EXAJVIPLES: 
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sage: find_rank_spread (2 , Graph ( { : [ 1 , 2 , 3 ] , 1 : [ 2 , 3 ] , 2 : [ 3 ] } ) ) 
(0, 1) 

780 sage: g = Graph ( { : [ 1 , 2 , 4 , 6 , 7] , 1 : [3 , 5 , 6 , 7 , 8] , ^ 

^2:[4 ,6 ,8] ,3:[4 ,7 ,6] ,4:[6] ,5:[6 ,8 ,7]}) 
sage: find_rank_spread (2 , g) 
( False , False ) 

V V 

subgraph=graph . copy () 
785 subgraph . delete .vertex ( vertex ) 

graph_bounds=minrank_bounds ( graph) 
if graph_bounds [0] = = graph_bounds [ 1 ] : 

# We have an actual min rank for graph 
subgraph_bounds=minrank_bounds (subgraph ) 
790 if subgraph_bounds [O]==subgraph_bounds [ 1 ] : 

# We have an actual min rank for the subgraph 
return graph_bounds [0] — subgraph_bounds [0] , <— > 
^subgraph_bounds [0] 
return False , False 



795 

def cut_vertex_connected_graph_mr ( c_vertex , graph) : 

?! 55 5? 

Given a cut vertex and a graph , attempt to calculate the minimum 
rank of the graph by applying the cut vertex method to the graph 
800 and vertex . 

INPUT: 

c_vertex — the cut vertex 

graph — the graph in which the cut vertex is contained 

805 

OUTPUT: 

a list of length 2 with the minimum rank as all entries , if 
the minimum rank can be calculated in this way 

810 False if the minimum rank cannot be calculated in this way 

EXAMPLE: 

sage : cut .vertex .connected _graph_mr (0 , Graph-*-^ 
^({0:[1,2,3],2:[3]})) 

(2, 2) 

815 sage: cut _vcr tex_connected _graph_mr ( 2 , Graphs 
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^({0:[1,2,3],2:[3]})) 
Traceback (most recent call last) 



ValueError : Supplied vertex is not a cut vertex 

55 55 5? 

820 g=graph.copy() 

if g . is_connected is False ■.#this should never happen 

raise ValueError , "Graph is not connected" 
if c_vertex not in graph . vertices () r^again , should never happen 
raise ValueError, "Supplied vertex is not in the graph" 
825 g . delete_vertex ( c_vertex ) 

subgraphs=g . connected_components_subgraphs () 

if len ( subgraphs ) <= 1: # C-vertex is not a cut— vertex 

raise ValueError, "Supplied vertex is not a cut vertex" 

830 

index=0 
rank_spread=0 
subgraph_mr_sum=0 
for subgraph in subgraphs : 
835 subgr aph_with_v = graph . subgraph ( subgraph . vert ices () 

^c_vertex ] ) 
new_rank_spread , subgraph_mr = \ 

find_r ank_spread ( c_vertex , subgraph_with_v ) 
if new_rank_spread is False: 
return False 
840 else : 

rank_spread += new_r ank_spread 
subgraph_mr_sum += subgraph_mr 

rank_spread = min ( rank_spread , 2) 
845 return subgraph_mr_sum+rank_spread , subgraph_mr_sum+r ank_spread 



def minrank_bounds ( graph , all_bounds=False , tests =[' precomputcd ' , 

■border', 'zero forcing', 'not path', 'forbidden minrank 2', 'not 
^planar', 'not outer planar', 'clique cover', 'cut vertex', '<-^ 
^disconnected', 'diameter']): 

55 55 55 

850 Find lower and upper bounds for the minimum rank of a graph. If 

all_bounds is False , then only return the best lower and upper 
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bounds. If True, return two dictionaries giving all applicable 
lower bounds and upper bounds, respectively. 



INPUT: 

graph — the graph whose minimum rank is bounded 

all_bounds — if False , then only return the best lower and <- 
"^upper 

bounds. If True, return dictionaries giving all 

^applicable lower 
bounds and upper bounds . 

tests — a list of tests to get bounds. Possible values are 
' precomputed ' , 'order', 'zero forcing', 'not path', 
'no forbidden', 'not planar', 'not outer planar', 'clique 

^ cover ' , 
'cut vertex', 'disconnected' 

OUTPUT: 

the lower and upper bounds for the minimum rank , in that 
border 



EXAMPLES: 

sage : minrank_bounds ( graphs . CompleteGraph (3) ) 
(1, 1) 

sage: minrank_bounds ( graphs . CompleteGraph (3) , all_bounds=True 

({'precomputed': 1, 'rank': 0, 'zero forcing': 1}, 

{' clique cover ' : 1 , 

' no forbidden ' : 2 , 

' not path ' : 1 , 

'order ': 2 , 

' precomputed ' : 1 , 

'rank': 3}) 

sage: minrank_bounds ( graphs . PathGraph (4) , all_bounds=True) 
({'cut vertex (1) ': 3, 'precomputed': 3, 'rank': 0, 'zero 
^forcing ' : 3} , 

{' clique cover ' : 3 , 
' cut vertex ( 1 ) ' : 3 , 
' order ' : 3 , 
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' precomputed ' : 3, 
'rank ' : 4 , 

'zero forcing (tree) ': 3}) 

sage: minrank_bounds ( graphs . PathGraph (4) , all_bounds=True , ^ 

^tests =['order ' , 'zero forcing ']) 
({'rank': 0, 'zero forcing': 3}, 

{'order ': 3, 'rank ': 4, 'zero forcing (tree) ': 3}) 
sage : minrank_bounds ( graphs . HeawoodGraph ( ) ) 
(8, 10) 

sage: minrank_bounds ( graphs . HeawoodGraph () , all_bounds=True) 

({'rank': 0, 'zero forcing': 8}, 

{ ' clique cover ' : 21 , 

'not outer planar ': 11, 

'not path ': 12, 

' not planar ' : 10 , 

' order ' : 13 , 

'rank': 14}) 

sage : minrank_bounds ( graphs . PetersenGraph () ) 

(5, 6) 

sage: minrank_bounds ( graphs . PetersenGraph () , all_bounds=True) 

({'rank': 0, 'zero forcing': 5}, 

{' clique cover ' : 15 , 

'not outer planar ': 7, 

' not path ' : 8 , 

' not planar ' : 6 , 

'order ': 9, 

'rank': 10}) 

V 55 ?5 

if isinst ance ( tests , str): 
tests = [ tests ] 

possible_tests = set ([' precomputed ' , 'order', 'zero forcing', 

^not path', 'forbidden minrank 2', 'not planar', 'not outer <— 
^planar', 'clique cover', 'cut vertex', 'disconnected', '-(-^ 
^diameter ' ] ) 

# Check tests 

unknown_tests = set ( tests ). difference ( possible_tests ) 
if len ( unknown_tests ) >0: 

print "Unknown tests specified: ", 1 i s t ( unknown_tests ) 

g=graph.copy 
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lower_bound = {'rank': 0} 
925 upper_bound = {'rank': g. order ()} 



if g . is_connected ( ) : 

bounds = min_rank_by_bounds ( graph , tests=tests) 
lower _bound . update ( bounds [ ] ) 
930 upper_bound . update ( bounds [ 1 ] ) 

# Try finding a cut vertex 
if 'cut vertex' in tests: 

c_vertex=find_cut_vcrtcx (g) 
936 if c_vertex is not False : 

cut_vertex_bounds = cut_vertex_connected_graph_mr 

^c_vertex , graph) 
if cut_vertex_bounds is not False : 

lower_bound [ ' cut vertex (%s ) '%(c_vertex , ) ] = 

^cut_vertex_bounds [0] 
upper_bound [ ' cut vertex (%s ) '%(c_vertex , ) ] = <^ 
^cut_vertex_bounds [1] 

940 else : 

if 'disconnected' in tests: 

connected-components = g . connected_components_subgraphs () 
lower_bound[ 'disconnected'] =0 
upper_bound [' disconnected ' ] = 
946 for component in g . connected_components_subgraphs ( ) : 

sub_bound = minrank_bounds ( component , tests=tests) 
lower_bound [' disconnected ' ] += sub_bound[0] 
upper_bound [' disconnected ' ] += sub_bound[l] 

960 # Make sure that the lower hound is not greater than the upper ^ 

'abound 

if max( lower_bound . values ) > min ( upper_bound . values ()) : 
raise StandardError , """ 
Best lower bound is greater than best upper bound; something is wrong-t 

lower bounds : %s 
966 upper bounds: %s" " " %(lower_bound , upper_bound ) 

if all_bounds is True: 

return lower_bound , upper_bound 
else : 
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960 



# Return the best lower and upper bounds 

return max( lower_bound . values ()) , min ( upper_bound . values () ) 



966 def wr ite_spreadsheet ( graph_numbers , filename): 

V 

Write a spreadsheet of minimum ranks and bounds for Atlas of 
■^Graphs graphs . 

DNfPUT: 

970 graph_numbers — a list of Atlas of Graph numbers 

filename — The filename for the spreadsheet. 
OUTPUT: 

975 the file is written (or overwritten if it already exists) 

EXAAIPLES: 

To write a spreadsheet of graph minimum ranks for all graphs 
in the Atlas of Graphs similar to the one that comes with 
^this 

980 program, first comment out the ' precomputed ' test in the 

min_rank_by_bounds function , then do 

sage: write_spreadsheet ( [ 1 . . 1 2 5 2 ] , " minrank . csv" ) 

7? V 

985 import CSV 

f = f i 1 e ( filename , "wb") 

fieldnames = ('Atlas number ',' Order ' , ' Size ', 'Minimum rank by ^ 
^program ' , 

'Program lower bound ',' Program upper bound', '<-^ 
^Connected ' , 

'Zero forcing LB', 'Diameter LB', 'Clique cover UB'< 
990 'Not planar UB','Not outer planar UB','Not path UB'< 

'Induced subgraph (minrank 2)', 'Cut vertex ',' Tree '< 

writer = csv . Diet Writer ( f , fieldnames=fieldnames ) 
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headers = {} 
995 for n in fieldnames : 

headers [n] = n 
writer . writerow (headers ) 



for graph_number in graph_numbers : 
1000 g = at las .graphs [ graph_number ]. copy ( ) 

lower , upper = minrank_bounds (g , all_bounds=True) 
row = {} 

row[ 'Atlas number'] = int (graph_number) 
row [' Order ' ] = g.order() 
1005 row ['Size'] = g.size() 



if max( lower . values ) = min ( upper . values ()) : 

row [' Minimum rank by program'] = max( lower . values () ) 
else : 

1010 row ['Minimum rank by program'] = None 



row 
row 



row 



row 
row 



row 
row 
row 



'Program lower bound 
Program upper bound 



max(lowcr . values () ) 
min ( upper . values () ) 



'Connected'] = g . is_connected () 

'Zero forcing LB'] = lower . get (' zero forcing', None) 
'Diameter LB'] = lower . get (' diameter ' , None) 



'Clique cover UB'] = upper . get (' clique cover', None) 
'Not planar UB ' ] = upper . get (' not planar'. None) 
'Not outer planar UB'] = upper . get (' not outer planar', 
^None) 

row ['Not path UB'] = upper . get (' not path', None) 



1025 # // we have a forbidden minrank 2 lower bound, then we have <^ 

^a 

# forbidden subgraph 

if 'forbidden minrank 2' in lower: 

row [ 'Induced subgraph (minrank 2) '] = True 
elif 'forbidden minrank 2' in upper: 
1030 row ['Induced subgraph (minrank 2)'] = False 

else : 

row [ 'Induced subgraph (minrank 2) ' ] = None 
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row [' Cut vertex'] = False 
for key in lower . keys () : 

if key . startswith ( ' cut vertex'): 

row [' Cut vertex'] = True 

break 

if 'zero forcing (tree) ' in upper: 

row [ ' Tree ' ] = True 
else : 

row [' Tree ' ] = False 
writer . write row (row) 
f . close ( ) 
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